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FINAL  REPORT 


GRANT  #:  NOQ0X4-O3-1-O212 

PRINCIPAL  INVESTIGATORS:  Drs.  Gareth  L.  Lawson,  Timothy  K,  Stanton  and 

Peter  H,  Wiebe 

INSTITUTION :  Woods  Hole  Oceanographic  Institution 

GRANT  TITLE:  Characterizing  Variability  in  the  Distribution  of  High- 
Frequency  Acoustic  Backscattering  in  a  Shallow  Water  Coastal  Region 

AWARD  PERIOD:  1  December  2002  -  31  August  2007 

OBJECTIVES;  The  central  purpose  of  this  GNR  Graduate  Traineeship  Award 
in  Ocean  Acoustics  was  for  the  Primary  Investigator,  Gareth  Lawson,  to 
design,  execute,  and  defend  his  doctoral  dissertation  research,  with  the 
Co-Primary  Investigators,  Timothy  K.  Stanton  and  Peter  H ,  Wiebe,  acting 
as  advisors.  In  terms  of  its  scientific  goal,  the  project  examined  the 
physical  and  biological  factors  underlying  variability  in  the 
distribution  of  high-frequency  acoustic  volume  backscattering  stemming 
from  zooplankton.  The  specific  research  objectives  were: 

1.  To  quantify  spatial  and  temporal  variability  in  zooplankton 
distribution,  and  hence  the  distribution  of  uncertainty  in  the 
acoustic  field,  in  shallow  water  coastal  regions, 

2.  To  assess  the  predictability  and  persistence  of  such  patchiness, 
and  understand  its  association  with  physical  and  biological 
oceanographic  processes . 

3.  To  continue  the  process  of  field-testing  and  refining  zooplankton 
scattering  models, 

APPROACH ;  The  overall  approach  was  to  apply  the  models  of  zooplankton 
scattering  developed  previously  at  Woods  Hole  Oceanographic  Institution 
( WHO I )  to  multi- frequency  volume  backscattering  data  collected  in 
continental  shelf  regions  under  other  research  projects,  in  order  to 
identify  dominant  scatterers  and  estimate  zooplankton  size,  abundance, 
and  biomass.  These  estimates  were  then  compared  between  times  and 
regions  to  assess  patterns  in  variability.  The  areas  under  study 
included  a  continental  shelf  region  west  of  the  Antarctic  Peninsula  and 
the  Gulf  of  Maine,  At  both  sites,  the  data  used  for  analysis  were 
collected  with  the  WHOI  Bio-Optical  Multi- frequency  Acoustical  and 
Physical  Environment  Recorder  (BIOMAFER-II ) . 

The  project  ultimately  came  to  focus  on  variability  in  the  distribution 
of  krill  (primarily  Euphausia  super ha  and  E.  crystallorophias)  in  the 
Antarctic  study  region.  An  important  aspect  of  the  research  involved  the 
development  and  testing  of  protocols  for  distinguishing  the  scattering 
of  krill  from  other  zooplankton  taxa,  delineating  krill  aggregations 
within  the  acoustic  record,  and  estimating  the  size,  abundance,  and 
biomass  of  aggregation  members.  These  methodologies  were  applied  to 
acoustic  data  collected  in  the  Antarctic  during  fall  and  winter  surveys 
conducted  by  the  Southern  Ocean  Global  Ecosystem  Dynamics  (GLOBEC) 
program.  The  resulting  estimates  of  krill  distribution  and  aggregation 
structure  were  then  related  to  features  of  the  physical  and  biological 
environment  measured  concurrently  during  the  surveys. 

ACCOMPLISHMENTS :  The  work  accomplished  under  this  award  made  important 
contributions  to  the  fields  of  both  zooplankton  acoustics  and  ecology, 
and  resulted  in  four  publications  with  the  Primary  Investigator  as 
senior  author  and  seven  as  co-author. 

Analyses  of  the  distribution  of  volume  backscattering  strength  at  120 
kHz  in  the  Antarctic  continental  shelf  study  region  revealed  strong 
seasonal  and  spatial  variability  during  the  falls  and  winters  of  2001 
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and  2002,  much  of  which  was  understandable  in  light  of  the  likely  effect 
of  meso-scale  circulation  on  the  distribution  of  zooplankton  scatterers 
(Lawson  et  al. ,  2004),  Four  general  types  of  scattering  features  were 
evident;  1)  large  bottom-associated  patches  found  in  coastal  regions  of 
variable  bathymetry,  2)  smaller  discrete  patches  found  primarily  in  the 
surface  mixed  layer  across  much  of  the  shelf,  3)  deep  diffuse  patches 
found  near  the  shelf-break,  and  4)  deep  homogeneous  scattering  layers 
situated  over  much  of  the  southern  portion  of  the  study  area. 

Predictions  of  expected  volume  backscattering  strength  based  on  depth- 
stratified  net  samples  and  taxon-specific  scattering  models,  in 
combination  with  analyses  of  Video  Plankton  Recorder  observations, 
suggested  that  the  first  two  feature  types  stemmed  from  krill  (Euphausia 
spp.)  aggregations  (Lawson  et  al,,  2004;  in  press  a).  Net  catches  and 
BIOMAPER-II  measurements  of  target  strength  suggested  that  myctophid 
fishes  were  the  dominant  scatterer  in  the  third  feature  type.  The 
available  evidence  concerning  the  deep  scattering  layers  was  less 
conclusive,  but  they  appeared  to  consist  of  a  complex  mixture  of 
zooplankton  scatterer  types,  perhaps  dominated  by  copepods  and  gas- 
bearing  siphonophores  (Lawson  et  al.,  2004), 

Similar  analyses  of  multi-frequency  volume  backscattering  data  (43,  120, 
200,  and  420  kHz)  collected  with  the  BIOMAPER-II  in  the  Gulf  of  Maine 
during  the  falls  of  1997-1999  under  projects  funded  in  part  by  QNR  and 
in  part  by  NOAA  revealed  strong  inter- annual  and  spatial  variability, 
both  of  which  differed  between  acoustic  frequencies.  Mean  backscattering 
levels  did  not  reflect  the  decrease  in  the  abundance  of  the 
ecologically-important  copepod  Caianus  finmarchicus  observed  in  net 
catches  in  the  fall  of  1998  relative  to  the  falls  of  1997  and  1999.  An 
increase  in  the  abundance  of  the  predators  and  competitors  of  C, 
finmarchiciis  may  have  served  to  keep  backscattering  levels  high  during 
this  period.  Distinct  patterns  and  diel  changes  were  evident  in  the 
vertical  distribution  of  backscattering,  although  these  differed  between 
basins  and  years.  Collaborative  work  with  Dr.  Andone  Lavery  of  WHOl 
suggested  that  at  certain  locations  in  this  Gulf  of  Maine  study  region, 
gas -bearing  siphonophores  dominated  scattering  at  all  four  BIOMAPER-II 
frequencies  (Lavery  et  al.,  in  press).  Otherwise,  no  single  zooplankton 
scatterer  type  dominated  volume  backscattering  at  all  frequencies 
considered. 

Following  these  initial  analyses,  the  present  project  came  to  focus  on 
krill  (primarily  Euphausia  superJba  and  E,  crystal lorophias)  at  the 
Antarctic  study  site.  This  decision  was  made  in  part  because  these 
animals  form  acoustically-identif iable  aggregations  within  which  they 
dominate  scattering  at  all  four  BXOMAPER-II  frequencies,  making  the 
estimation  of  biological  parameters  like  krill  abundance  and  size  from 
four -frequency  data  a  relatively  straightforward  problem.  The  decision 
also  related  to  the  fact  that  the  Antarctic  field  work  began  at  the  same 
time  as  the  Primary  Investigator  was  formalizing  his  dissertation 
research  plan.  Focusing  on  the  Antarctic  study  site  thus  allowed  him  to 
be  involved  in  both  data  collection  and  analysis,  unlike  the  Gulf  of 
Maine  study  site  where  the  data  were  collected  some  years  earlier. 

Much  of  the  initial  work  of  the  project  made  use  of  existing  scattering 
models  and  parameter! zat ions  for  interpreting  observed  backscattering 
levels.  In  the  case  of  Antarctic  krill,  however,  an  investigation  was 
made  into  improving  the  parameterization  of  a  Distorted- Wave  Born 
Approximation  (DWBA)  based  scattering  model  (Lawson  et  al,,  2006). 

Direct  measurements  made  with  a  Video  Plankton  Recorder  were  used  to 
constrain  the  key  parameter  governing  the  orientation  of  the  animal 
relative  to  the  incident  acoustic  wave.  Recently  published  length -based 
regressions  were  used  to  constrain  the  animal's  acoustic  material 
properties,  rather  than  the  earlier  approach  of  using  single  parameter 
values  for  all  lengths.  The  newly-parameterized  model  produced 
predictions  of  krill  target  strength  more  consistent  with  in  situ 


measurements  than  earlier  parameterizations .  These  predictions  were 
smaller,  however,  than  expected  under  the  semi -empirical  model 
traditionally  used  to  estimate  krill  target  strength*  Use  of  this  semi- 
empirical  model  in  estimating  krill  abundance  would  result  in  estimates 
too  low  by  a  factor  of  at  least  2*75. 

with  this  fully- parameter! zed  and  field-tested  scattering  model  in  hand, 
methods  were  next  developed  and  refined  for  identifying  krill 
aggregations  and  estimating  krill  length,  abundance,  and  biomass  on  the 
basis  of  multi -frequency  acoustic  data  alone,  without  recourse  to  net 
samples  (Lawson  et  al.,  in  press  a)*  A  biologically-based  threshold 
level  of  volume  backseat ter ing  strength  for  identifying  krill 
aggregations  was  derived  on  the  basis  of  krill  visual  acuity  and 
estimates  of  the  minimum  density  of  krill  over  which  a  given  individual 
could  maintain  a  visual  association  with  its  nearest  neighbor,  and 
thereby  with  the  aggregation  as  a  whole*  Differences  in  mean  volume 
backseat ter ing  strength  at  120  and  43  kHz  further  served  to  distinguish 
krill  from  other  sources  of  scattering*  A  multi-frequency  inverse  method 
was  developed  to  estimate  the  mean  length  and  numerical  density  of  krill 
in  these  acoustically- identified  aggregations  on  the  basis  of 
measurements  of  mean  volume  backscattering  strength  (43,  120,  200,  and 
420  kHz)  *  The  biomass  density  of  krill  was  then  estimated  based  on 
measurements  of  volume  backscattering  strength  at  120  kHz,  the  mean 
length  of  krill  estimated  from  the  multi -frequency  inversion,  the  target 
strength  in  relation  to  length  model  of  Lawson  et  al .  (2006),  and  the 

wet  weight -to- length  relationship  measured  by  Wiebe  et  al  *  (2004)* 

Measurements  were  also  made  of  certain  features  (e*g.,  position,  shape, 
internal  structure)  of  each  identified  krill  aggregation. 

These  methods  were  tested  at  locations  where  independent  information  on 
the  presence  of  krill  as  well  as  on  krill  length  and  abundance  were 
available  from  Video  Plankton  Recorder  observations  and  net  catches  from 
a  1-rrf  Multiple  Opening /Closing  Net  and  Environmental  Sensing  System 
(MOCNESS) *  Comparison  of  the  results  of  the  inverse  analysis  to  net 
samples  were  favorable  in  the  case  of  estimated  mean  krill  length 
(within  a  margin  of  error  of  3-12%).  Acoustic  estimates  of  krill 
numerical  density,  however,  exceeded  those  from  nets  by  one  to  two 
orders  of  magnitude,  likely  due  primarily  to  net  avoidance  and 
secondarily  to  differences  in  the  volumes  sampled  by  the  two  systems. 

The  potential  for  multi- frequency  data  and  mathematical  inverse 
techniques  to  be  used  for  the  simultaneous  and  quantitative  estimation 
of  zooplankton  abundance  and  size  has  been  known  for  some  time,  but  to 
the  best  of  our  knowledge,  this  project  marks  the  first  time  that  such 
methods  have  been  applied  to  broad- scale  data  from  Antarctic  krill 
surveys  * 

Application  of  these  methods  to  multi- frequency  survey  data  collected 
during  the  fall  and  winter  2001  and  2002  surveys  at  the  Antarctic  study 
site  with  the  BIOMAPER-II  demonstrated  strong  seasonal,  inter-annual, 
and  spatial  variability  in  the  distribution  of  krill  (Lawson  et  a!*,  in 
press  b) .  The  distribution  of  krill  aggregations  was  characterized  by 
many  small  aggregations  closely  spaced  relative  to  one  another, 
punctuated  by  much  fewer  aggregations  of  very  large  size  that  accounted 
for  the  majority  of  overall  biomass  in  the  region.  Generalized  additive 
models  indicated  that  high  krill  biomass  was  consistently  associated 
with  regions  close  to  land  where  the  temperature  maximum  below  200  m  in 
depth  was  cooler  than  was  available  over  the  shelf  as  a  whole*  Krill 
were  thus  not  associated  with  regions  where  intrusions  of  warm  and 
nutrient-rich  circumpolar  deep  water  onto  the  shelf  were  present, 
despite  such  intrusions  being  thought  to  be  an  important  driver  of 
primary  productivity  in  the  region. 

Finally,  the  morphology,  internal  structure,  and  vertical  position  of 
individual  krill  aggregations  were  examined  in  relation  to  a  variety  of 


properties  of  the  physical  and  biological  environment*  Krill 
aggregations  were  observed  to  exhibit  diel  changes  in  vertical  position 
and  biomass  density,  with  aggregations  deeper  and  denser  during  the  day 
than  during  the  night*  Such  diel  vertical  migrations  had  not  previously 
been  observed  for  krill  during  winter*  Concurrent  observations  of 
chlorophyll  a  concentrations  (an  index  of  phytoplankton  abundance)  and 
the  occurrence  of  predators  together  reinforced  the  conclusion  that 
aggregation  and  diel  vertical  migration  represent  strategies  to  avoid 
visual  predators,  while  also  allowing  the  krill  access  to  shallowly- 
distributed  food  resources. 

In  addition  to  the  main  dissertation  research  project,  the  Primary 
Investigator  was  involved  in  a  number  of  collaborations,  examining  the 
distribution  of  higher  predators  in  the  Antarctic  (including  seabirds, 
seals,  and  whales)  in  relation  to  their  krill  and  zooplankton  prey 
( Friedlaender  et  al*,  2006;  Ribic  et  al*,  in  press),  as  well  as 
examining  the  distribution  and  behavior  of  Antarctic  krill  and  other 
zooplankton  using  nets  (Ashjian  et  al*,  2004;  Wiebe  et  al.,  2004)  and  a 
Video  Plankton  Recorder  {Ashjian  et  al, ,  in  press). 

CONCLUSIONS :  1*  At  certain  acoustically-identif iable  times  and  locations 

in  the  Antarctic  study  region,  krill  were  the  dominant  scat ter er  at  all 
four  BXGMAPER-II  frequencies  (43,  120,  200,  and  420  kHz)*  This  makes  the 
krill  a  favorable  species  for  study  with  multi -frequency  acoustics. 

2*  Away  from  these  locations,  other  zooplankton  dominated  scattering, 
especially  copepods  and  gas-bearing  siphonophores *  This  suggests  that 
the  assumption  often  made  in  acoustic  surveys  for  Antarctic  krill  that 
all  scattering  stems  from  krill  is  inappropriate,  and  care  must  be  taken 
to  identify  properly  those  exact  locations  where  krill  are  responsible 
for  observed  scattering* 

3  *  Strong  spatial  and  inter-annual  variability  in  volume  backscattering 
was  observed  at  the  Gulf  of  Maine  study  site*  Ongoing  analyses  are 
addressing  the  causes  of  these  inter -annual  changes* 

4.  Antarctic  krill  were  observed  to  orient  themselves  mostly 
horizontally,  or  at  normal  acoustic  incidence  relative  to  a  vertically- 
oriented  echosounder*  Parameterization  of  a  DWBA-based  scattering  model 
for  krill  with  the  observed  distribution  of  angles  of  orientation  and 
direct  measurements  of  all  other  necessary  parameters  led  to  target 
strength  predictions  highly  similar  to  in  situ  measurements* 

5*  Multi- frequency  inversions  can  be  used  successfully  to  estimate 
Antarctic  krill  length  and  abundance  over  large  survey  areas  without 
recourse  to  net  samples, 

6.  Highest  Antarctic  krill  biomass  during  fall  of  both  survey  years  was 
consistently  associated  with  regions  close  to  land  where  temperatures 
below  200  m  in  depth  were  cooler  than  what  was  available  over  the  shelf 
as  a  whole. 

7*  Krill  aggregation  and  diel  vertical  migration  represent  strategies 
to  avoid  visual  predators  during  the  day,  while  also  allowing  the  krill 
access  to  shallowly-distributed  food  resources  by  night. 

SIGNIFICANCE :  This  project  has  resulted  in  a  detailed  understanding  of 
variability  in  the  distribution  of  zooplankton  volume  backscattering, 
particularly  in  the  case  of  Antarctic  krill.  Ultimately  this  should 
allow  such  variability  to  be  modeled  and  predicted*  This  constitutes  an 
essential  step  in  constraining  the  uncertainty  introduced  into  Navy 
representations  of  the  acoustic  field  by  these  important  scatterers,  and 
in  understanding  the  implications  of  zooplankton  patchiness  to  the 
distribution  of  higher  predators,  including  marine  mammals  and  exploited 
fish  species*  Much  of  the  present  research  has  made  use  of  existing 
scattering  models.  Through  comparisons  of  acoustic  data  to  video  and  net 
samples,  an  approximate  model  for  Antarctic  krill  has  been  further 
evaluated  and  better  parameterized.  The  outcome  is  a  field-validated 
acoustic  scattering  model  of  a  complex,  naturally  occurring  scatterer* 


PATENT  INFORMATION:  None. 
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